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Abstract 

We develop a supersymmetric representation of spin operators which unifies 
the Schwinger and Abrikosov representations of SU(N) spin operators, allow- 
ing a second-quantized treatment of representations of the SU(N) group with 
both symmetric and antisymmetric character. By applying this to the SU(N) 
Kondo model, we show that it is possible to develop a controlled treatment 
of both Magnetism and the Kondo effect within a single large N expansion. 
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I. MOTIVATION FOR A NEW SPIN REPRESENTATION 



Recent experiments on quantum phase transitions in heavy fermion materials have led 
to a debate about how magnetism condenses out of the metallic state at absolute zero. Cer- 
tain heavy fermion materials can can be tuned between the magnetic and the paramagnetic 
state through the use of pressureBJl or chemical dopingi'S. The quantum critical point which 
separates these two phases is of great current interest, in part because materials in its vicin- 
ity may become fundamentally new kinds of metaH Heavy fermion materials contain a 
dense lattice of magnetic moments; conventional wisdom assumes that the spins of the local 
moments are magnetically screened and of no importance to the magnetic quantum critical 
point!~0. Recent neutron data contradict this viewpoint, by showing that the spin correla- 
tions at the quantum critical point are critical in time, but local on an atomic scaleBElHll, 
suggesting that unscreened local moments emerge from the metallic state at the quantum 
critical point. 

If it is indeed true that the magnetic quantum critical points involve local moment 
physics, then a new theoretical approach is required. Traditionally, heavy fermion physics is 
modeled using a Kondo lattice Hamiltonian,0 describing the interaction between a bath of 
conduction electrons and an array of local moments. One of the well-developed theoretical 
methods for approaching this model is the large N expansionllHi!, where the idea is to use 
a generalization of the quantum mechanical spin operators, in which the underlying spin 
rotation group is generalized from SU(2) to SU(N). The utility of this method derives from 
the fact that in the limit iV — > oo, it provides an essentially exact, analytic treatment of 
the Kondo lattice problem. Unfortunately the way this procedure is carried out at present, 
magnetic interactions are suppressed as a 1/N 2 correction, beyond the horizon for a con- 
trolled computation. In this paper, we show how we can overcome this shortcoming by the 
use of a supersymmetric spin representation for local moments. 

The theoretical description of interacting local moments poses a fundamental problem: 
the Pauli Spin operator S does not satisfy a Wick's decomposition theorem, which pre-empts 
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its use in a Feynman diagram approach. The traditional solution is to represent the spin 
in terms of either bosons, or fermions. In the "Schwinger boson" approachJUill the spin 
operator is represented in terms of an iV-component boson b a ; in the alternative "Abrikosov 
pseudo-fermion" representation,!^! the spin is represented by an N component fermion, f a , 
as follows 



(a) Antisymmetric 



(b) Symmetric 



Q 



nf — Q 



2S 

n b = 2S 



Fig. 1. Young tableauxEl for (a) antisymmetric and (b) symmetric representations generated 
by the Abrikosov fermion and Schwinger boson representations respectively. 

Here, T = (T 1 , . . .T M ) represents the M = (N 2 — l)/2 independent SU(N) generators. By 
combining 2S "Schwinger bosons" together, one generates a symmetric representation of 
SU(N), denoted by a horizontal Young- Tableau with 2S boxes (Fig. lb). Conversely, in the 
pseudo-fermion approach. Q < N spin fermions are combined to generate an antisymmetric 
representation of SU(N), denoted by a column Young- Tableau with Q boxes, (Fig. la). 

Most cases of physical interest correspond to an elementary spin, or a single box in the 
Young- Tableau. Unfortunately, to develop a controlled Feynman diagram expansion, we are 
obliged to consider a large number of such boxes, letting iV — > oo keeping either q = Q/N or 
m = 2S/N fixed. In the process of letting iV — > oo, some essential physics is lost. Symmetric 
representations are ideal for treating magnetism, where the ordered moment involves a highly 
symmetric condensation of spin bosons, but they lose all information about the Fermi liquid 



fixed point. Antisymmetric representations capture the development of the Kondo effect 
and heavy fermion bands, but magnetism is suppressed. 

Various authors have tried to develop alternative representations of the spin operatorBil 
and in this context one interesting idea is to use supersymmetry to simultaneously express 
the bosonic and fermionic character of local moments0~il. 



II. SUPER SPINS 

We now examine a class of spin representations which preserve both symmetric and 
antisymmetric correlations. Consider the spin operator that is a sum of n/ fermions and rif, 
bosons, given by 

S = S / + S 6 (1) 

By combining Q = rif + bosons and fermions together, we generate "L-shaped" represen- 
tations of SU(N). For each choice of rif and n?,, we generate two irreducible representations. 
For example, we can combine one fermion and two bosons as follows: 



f 




b 


b 



b 


b 


f 





"S = 1/2 



l S = 3/2" 



For SU(2), these correspond to a spin 1/2, and a spin 3/2 representation. To uniquely pa- 
rameterize an irreducible representation, we need to fix the Cazimir S 2 = J2 a S a S a - Consider 
an L-shaped representation of SU(N) of width w, height h : 



' w 
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If the generators of the fundamental representation are normalized according to Tr|F a r 6 ] = 
5 ab , then the expression for the Cazimir of an arbitrary irreducible representation isllll 

s 2 = Q(N ) l ~ Q) + e ^k+i- 2j ) ( 2 ) 

JV j=l,h 

where rrij is the number of boxes in the j — th row from the top and Q is the total number 
of boxes. For an L-shaped Young tableau, (mi,m 2 • • • m h) = (w, 1, 1 . . . , 1), so that 

S 2 = g(iV ^~ Q) + w{w - 1) - h{h - 1) (3) 

If we substitute Q = w + h — 1 and Y = h — w we then obtain 

S 2 = Q(N -Y -Q/N). (4) 

In this way, each irreducible L-shaped representation of SU(N) is uniquely defined by the 
two quantities (Q,Y), where Y can assume the values 

Y = -Q + 1,-Q + 3,...Q-1. (5) 

For example, if Q — 3, there are three irreducible representations: 

Y= 2 -2 



We now seek to cast both Q and Y in an operator language. In terms of the boson and 
fermion operators, the Cazimir can be written 

S 2 = (S / + S,) 2 (6) 

If we expand this expression (appendix A) by using the completeness relation 

-Wr^r^ = N5 ai 5fS8 - 8 a p8 Sl (7) 



we are able to express the Cazimir in operator form: 



S 2 = Q(N-y 



Q 



(8) 



N 



where now 



Q = n f + n b , 



(9) 



fixes the number of boxes and 



v 



y = nf - nb+ [0,0t] 




(10) 



is the operator measuring the asymmetry h — w of the representation. Here we have intro- 
duced the operators 



If we wish to study a spin system described by the (Q, Y) representation, then we must 
restrict our attention to states in the Hilbert space which satisfy 



Curiously, although this constrains the total number of bosons and fermions, the difference 
n^ — nf is only partially constrained, reflecting the fact that bosons and fermions can inter- 
convert without altering the representation. 

When we represent spins in terms of bosons and fermions, each box in the Young tableau 
is associated with a fermion or boson, where fermions occupy the column, bosons the row. 
The corner of the tableau can contain either a boson, or a fermion. The operator 6^ converts 
the corner box from a boson into a fermion, whilst 9 converts it back again. These operators 
are the generators of a "supergroup" S77(l|l)0, with the algebra 



= blfa- 



(11) 




(12) 



{9\9} = Q. 



(13) 
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The spin operator commutes with these generators 

[S 6 + S f ,0] = [S 6 + S,,0t] = O, (14) 

so that the representation is supersymmetric. Q and 3^ also commute with 6 and 6\ and 
these are the Cazimirs of this group. From QT3|), we see that operators P b = -^09^ and 
pf = ^0^0 satisfy Pb + Pf = 1: they are the projection operators which respectively project 
out states with "bosons" or "fermions" in the corner of the Young tableau. In this way, we 
see that V = Pf — Pb is +1 or —1, depending on whether the state has a boson, or fermion 
in the corner of the representation For example, the representation given by (Q, Y) = (3, 0) 
can be written in two ways: 



b 


b 




f 


b 


f 
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f 





V = 1 V = -1 
n/ — rib — — 1 rif — rib = I 

Y = Y = 



The invariance of the representation under the boson-fermion transformation is a manifesta- 
tion of the supersymmetry. When a boson is converted into a fermion, the change in rif — rib 
is compensated by the change in V, so that y is invariant. 

We can alternatively write the constraints in terms of the height h = (Q + y + l)/2 and 
width w = (Q - y+ l)/2 of the tableau: 



2S 




For the fundamental representation, described by the state 
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n*f — 2S — 1. Independently of the way we represent the spin, some bosonic and fermionic 
character is always present, reflecting the fact that a spin can give rise to a "bosonic" local 
moment, or it can produce "fermionic" singlet bound-states. Traditionally, one of the above 
constraints is dropped: in the approach now adopted, both constraints are simultaneously 
applied. 

There are two ways in which we can use the new constraint. We can work within a 
"grand canonical" ensemble, where Q is fixed, but j) is associated with a chemical potential, 

H' = H + (y (17) 

By tuning ( from negative, to positive values, the ensemble is driven from an antisymmetric 
to a symmetric representation. In fact, since the cazimir S 2 = Q(N — y — ^) is linearly 
related to y, a finite value of ( is physically equivalent to the introduction of a Hund's 
interaction into the Hamiltonian. 

H' = H-^(S 2 , (18) 

where a constant term ((N — Q / N) has been omitted. The supersymmetric spin represen- 
tation thus enables us to progressively increase the strength of the magnetic interactions by 
tuning the spin representation. 

Alternatively, we may work with a definite representation, where y — Y . The partition 
function for this model is 

Z[Q ,Y} = Tr[P QoX e^ H } (19) 

where Pq D) y projects out the states with definite Q = Q Q and y — Y. By specifying 
these two constraints, we are still working in an ensemble where the individual number of 
fermions or bosons are not separately constrained, and in this way, we are able to develop 
a supersymmetric field theory. We can implement these two constraints by carrying out a 
Fourier transform over the chemical potentials A and ( associated with Q and y respectively, 
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Z[g °' Y]= I -^f^[e- p{H+X ^- QMS - Y) ] (20) 
where both ( = X a +ix and A = \ +iy are integrated along an imaginary axis, x, y G [0, 2wT]. 



III. APPLICATION TO THE UNDER-SCREENED KONDO MODEL 

A. Formulation of Lagrangian 

To illustrate the approach, we develop it for the single impurity Kondo model, given by 

Hp Hk Hq^ ^-S—, 

h = E *h<L<** + jr<iF°fiC0 -S + XQ + j-Qy (21) 

k,a 

Here, H describes the conduction electron sea, Hk is the interaction between the conduction 
electron spin density, and the local moment, where = n s -1 / 2 J2k c ka creates an electron 
at the site of the local moment (n s = no. of sites). Hq and Hy impose the constraints. 
(Note the way in which Hy has been written: by multiplying the operator y by Q, we cast 
it in a form which is unchanged upon normal ordering: Qy =: Qy :. By writing it in this 
way we can immediately translate the fields to their coherent state representation. Inside a 
path- integral, where Q = Q is imposed, we are then able to replace rrQy —> y.) 

To date, the one-channel Kondo model has been studied in the large N approach for 
completely symmetricil and completely antisymmetric^ representations of S. For the latter 
the local moment is quenched to form a local Fermi liquid; for symmetric representations 
the spin S is only partially screened by the Kondo effect to form a spin S — ~. By tuning 
the representation from the one limit to the other, we are able to examine how the local 
Fermi liquid interacts with the emergent local moment as the local moment grows. For 
fully symmetric representations, it is known that the residual Fermi liquid decouples from 
the partially screened momentiil. One of the surprising discoveries of this study, is that in 
intermediate representations, the heavy Fermi liquid and the partially screened moment can 
become antiferromagnetically coupled and the strong-coupling fixed point becomes unstable. 
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We shall see that in this simple model, a two-stage Kondo effect then takes place; richer 
consequences are likely in a lattice model. 

Our first step is to write the constrained partition function as a path integral 

Z = J V[c, f, b, A, Q e -Io[ C o+H K +C 3Usy -(XQ +CY)}dr ^ 

where we have divided the action into three terms: 

Co = 5I c L(<9t + e k )c kCT , 

Csusy = Y\fi d rU + + H Q + Hy. (23) 

The first term describes the conduction electrons, in the second term we have rewritten 
the local spin in terms of slave-fields, and the third term contains the machinery of the 
super symmetric representation. We use a single notation for the field operators and the 
c- number fields that represent them inside the path integral. 

Our next step is to formulate the Lagrangian in a form that clearly exhibits the su- 
persymmetry. We shall begin by casting C susy in a form which is gauge invariant under 
time-dependent super-rotations. It is convenient to combine the slave fields into a single 
spinor, 



Using this notation, 



Csusy = E H 9 T + A + (T 3 V a - ^6 (25) 



where T3 is a Pauli matrix. Since the starting Hamiltonian and each of the constraints 
commutes with the super-generators, the full Lagrangian is invariant (appendix B) under 
time-independent super-rotations ^ a — > g^ a where 
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(26) 



-fj V 1 - wj. 

is an element of the supergroup SU(1\1). If we make this transformation time dependent, 
the derivative terms become 



(27) 



Expanding the second term, we obtain 



tf a {gtd T g)* ff = tfd rV + fjd T e + Q fjd TV 



(28) 



where we have replaced ^ a ^ a — > Q inside the path integral. Since Z is unchanged by this 
change of basis, we can integrate over all g{r) 



Z = J D[fj,rj\ J D[c, f i b ) \X]e-^ c+e ' d ^ dTe+Q ° nd ^ dT 
= jD[c, f, b, A, C] expj- J\c - ^-fd T 6]dr] 



(29) 



By absorbing the additional term into a redefined 



susy 



d T + A + Cr 3 



m a - e\d T + 2Q6 

Wo 



(30) 



the Lagrangian becomes invariant under time-dependent super-rotations. The first term in 
C* susy describes the level splitting between the bosonic and fermionic components of the spin. 
The second term describes a residual interaction between the spin and heavy electron fluid. 
We can factorize this term, to obtain 

Hi 

-6\d T + 2Q0 - Q Q a\d T + 2Qa + [f%(d T + 2()a + a\d T + X^U] (31) 
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The first term tells us that the field a represents a dynamical fermion with the commutation 
algebra {a, a'} = 1/Q ■ This spin-less particle mediates the interaction between the spin 
and the Fermi liquid; Hj defines the vertex for the decay process /~ # b a + a~ . To represent 
this process, we denote the propagator for the a fermion by the Feynman diagram 



n 



(32) 



The vertices which inter-convert the heavy electron and spin bosons will be denoted by 



> =i(d T + 2() = i(2(-iu; n ), 



(33) 



where the "i" is required to give the correct amplitude for the exchange of the gauge fermion 
and 



G b (iis n ) = {iv n - A 6 ) 1 
Gf(iu n ) = {iu n - Xf Y 1 



(34) 



represent the propagators for spin bosons and the f-electrons. The mediated bare interaction 
between the spins and the heavy f-electrons is then 




r a Vo 



(35) 



It is a rather unique feature of this kind of approach that the spin interactions are carried 
by fermions, rather than bosons. Our final form for C* usy can now be compactly written 

where we have defined (d T + 2() = 0, (—d T + 2() = 0. 

As our next step, we factorize the Kondo interaction term Hk 



IT ^ 



^ a + Q o a*0a 



(36) 



ctVf a + Vflc a ) + (01^ + ^) 



(37) 



where V is a complex c-number field and is its fermionic partner. If we now introduce 

V 



V 



(38) 



then the transformed Lagrangian takes the form C = C Q + C* + H* K) where 
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£ Q = $3 c L(<9t + e k )c kCT , 

ko- 



iV. 



IT " 



Let us briefly examine the gauge invariance of this Lagrangian. If 



(39) 



(40) 



is a general S77(l|l) matrix, where g takes the form (^), then under the gauge transforma- 
tion — > /i^o-, V — > 7iV„ £ Q and ii# are invariant, but becomes 

~ C 
)at -C 

When we expand the first term (Appendix B), we find that 



h^ a + Q Q a ] $a 



(41) 



(42) 



where A' = A + C' = C + anc ^ = + ?y)e 2i6 * c , so C is invariant under the gauge 
transformation, 



A -> A - z0 Q) C - C - % 



a — > e~" c a — 7], a' — > e """'•a 1 — 77. 



(43) 



This gauge invariance leads to bosonic and fermionic zero-modes. To eliminate them, we 
must carry out a gauge fixing procedure. We can always parameterize V in the form 





~v ~ 






V = h 




= 9 















(44) 



or written out explicitly, 



V 



V n e lt) f 



-fj 



(45) 
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where Of = 0q + 9$ and V = \JVV + cjxj) is real. This transformation uniquely specifies 
both fj = —y4> an d the phase factor and e ldf = ■Xr, but does not specify 9 b = 9q — 9^. We 
shall adopt a gauge choice where 9b = 0. By applying the gauge transformation (f|^), we 
can absorb the fermionic fluctuations in V into a redefinition of the fields. The gauge fixed 
hybridization is now 

r 1 NV 2 

H* k = V Y: clf. + flc +—^. (46) 

a 

With our gauge choice 9 b = 0, the variable A/(r) = A + ( + i0f becomes dynamical , but the 
variable A& = A — ^ + i9b = A — C remains a time-independent integration variable. In this 
gauge fixed form of the Hamiltonian, the interaction between the Fermi and spin fluids is 
entirely contained within £* usy , and it is here where we should look if we are to obtain new 
physics. 

In the existing large N approach to the Kondo model, magnetic interactions are a 1/N 2 
correction to the mean-field theory, an order of accuracy that is beyond current theoretical 
approaches. A supersymmetric approach enhances the magnetic interactions by a factor of 
iV 2 , bringing them within the realm of Gaussian fluctuations about a new large N mean- 
field theory. To carry out concrete calculations, we expand around a large iV saddle point 
of the path integral obtained by taking N oo, maintaining Q/N = q and Y/N fixed. By 
allowing the number of bosons n b = Nn^ to become large, the Bose field is able to condense 
and form a magnetic moment: 



(K) = ^Mz a ~ 0(VN) (47) 

where z a is a unit spinor. The magnetic interaction between spins at different sites is given 
by RKKY diagrams (Fig. 0). The factors of y/N associated with the Bose condensate 
produce an the inter-site magnetic interaction of order 0(1) : a factor of iV 2 enhancement. 
These magnetic corrections appear as part of the Gaussian fluctuations of the a fermion and 
by calculating them, we are able to carry out a controlled treatment of magnetism and the 
Kondo effect. 
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FIG. 2. Magnetic interaction between two spins at site i and j within the current approach. 
(Note that when V a > 0, the f-fermion can propagate from site to site.) The factor (VW) 4 associated 
with the Bose condensate, enhances the magnetic interaction by a factor of N 2 , so that it appears 
as the as the first in a series of RPA diagrams associated with the Gaussian fluctuations of the a 
fermion. 

B. Mean Field Theory 

To illustrate this kind of calculation, we develop the machinery for the single impurity 
model. Although there are no inter-site magnetic interactions, the machinery of the su- 
persymmetry is needed to compute the magnetic interaction between the partially screened 
local moment and the Fermi liquid in the single impurity model. The techniques that we 
now illustrate can be generalized to the lattice. 

Our first step is to formally integrate out the conduction(c) and the slave fields (/, b) 

Z = J D[a,V,\,(]e- S °" 
e -s eff = J D[c, f, b]e~ Jo( c °+ c ™ S y+HK)dr ( 4g ) 

Since the second integral is bilinear in the fields, it can be carried out to yield 
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STr 



Tr[F] - Tr[B] (50) 



S eff [a, A, C] = -STrln[<9 r + A + (r 3 + E] + f \Q a^a - \Q - (Y]dr (49) 

Jo 

and S is the self-energy correction induced by the coupling to the conduction electrons and 
the a fermion and STr[. . .] denotes the "super trace"- the difference between the fermionic 
and the bosonic trace: 

" F a" 
J B_ 

Our procedure is then to expand the effective action to quadratic order around the saddle 
point where A, ( and V(r) = V are static, and a — 0. 

S eff [V, A, C, a] = S MF + 0(5 A 2 ) (51) 

where S A denotes the fluctuations. S M f = (3F MF determines the leading O(N) mean- field 
contribution to the Free energy. By carrying out the Gaussian integral over the fluctuations 
in SA we can then determine the 0(1) correction to the mean-field theory and the magnetic 
interactions within the medium. 

We begin by computing the saddle point, described by the the mean-field Hamiltonian 

H mft = ]T e k cLc kCT + Vo^Hfa + ftc a ] + (A - CK + (A + C)n f - \Q a (52) 

ko- <r 

where V Q and A are to be determined self-consistently. This mean-field theory describes a 
Kondo resonance formed between the conduction electrons and the antisymmetric part of 
the spin. The residual symmetric part of the spin is unquenched, and described by a sharp 
bosonic level at energy A^ = A — (. For the moment, we shall work in the ensemble of definite 
(, examining how the mean-field evolves as we increase ( to favor representations that are 
increasingly symmetric. If we define A/ = A + (, A& = A — (, then in the presence of the 
finite hybridization, the Greens functions for the slave fields are now given by 



= G b {iu n ) = (%v n - A 6 ) 1 

= G f (iu n ) = (iu n - \ f + iAsign(cu„)) _1 (53) 
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where A = nV^p is the hybridization width of the conduction electron and p is the conduc- 
tion electron density of states. The mean-field Free energy is then given by 

F mft = NTY,{\n[-Gl\iu n )] - ln[-Gj 1 (ic < ; n )])e <w » 0+ + N^f - XNq (54) 

The first term is just the free energy of a free boson. Carrying out the frequency sum on 
the second term, (Appendix C) we obtain 

%^ = + Mh) + V ~f - M (55) 

where 



$ f ( z ) = -2TReln 



F /Z±L 



D/2 



<S> b {z)=T\n{l-e- pz ) (56) 

are the fermionic and bosonic contributions to the energy, D is the conduction band half- 
width and £ = A/ + iA. If we differentiate this result with respect to A and A, we obtain 

7i fdF .dF\ / f 1\ _ / £> \ 1 
N{dA +t ^x) =^{^ + 2)- ln {2?nT) + Tp +miq - nb) 

■^+ 5 )-K JS Sff-)=° (57) 



2thT 2 

where ^(z) = d z \nT(z) is the digamma function, n fe = ra^/iV = [exp(Afe/T) — and 
T K = _Dexp[— 1/Jpo] is the Kondo temperature. At zero temperature we may replace, 
■ip(z) — > Inz, so the T = mean-field equations are then 

£ = 7^(9-**) (58) 

If the Bose field condenses, then A& = 0, so A = (, A/ = 2( and £ = 2( + iA. In this case, 
we can solve for the size of the unquenched moment, M = Nm = Nh^ and the width A of 
the Abrikosov Suhl resonance with which it co-exists: 

m = q-- cos' 1 — , (C > Cc) 



7T V T) 



A = ^(7^ - (2C) 2 , (59) 
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FIG. 3. Top: magnetization of partially screened local moment as a function of £. Bottom: 
Phase diagram for the supersymmetric impurity Kondo model, showing how the representation 
of the local moment evolves as ( is increased. Shaded area indicates under-screened region. For 
C < 0, the under-screened region involves a two-stage Kondo effect. 

where 

C c = (7V2)cos(7rg). (60) 

corresponds to a critical value of the Hund's interaction beyond which the local moment 
develops an unscreened component. There are three regions: 

q 2(/T K > 1 

? " \ 008-1 (ft ) > 1 > 2 C/7V > cos(Trg) (61) 
cos(Trg) > 2(/T K 

corresponding to an unscreened, partially screened and fully screened local moment. Fig. 1. 
shows the mean-field phase diagram. 

Next, let us consider the residual interactions between the un quenched local moment, 
and the Fermi liquid. The constraint term (y generates the residual interaction 



m 



Hi = -f{f%){blU>) (62) 

between the heavy-f electron and the unscreened moment. Conventional wisdom, based on 
the spin S Kondo model supposes that such residual interactions are always Ferromagnetic. 
For C > 0, this is indeed the case, and the fixed point described by the mean-field theory is 
thus stable. By contrast, if £ < 0, then the residual interaction is antiferromagnetic. The 
presence of such terms is unexpected. We shall see that for the impurity, this leads to a 
two stage Kondo effect. In the single impurity model, the condition ( < corresponds to 
the requirement that n*j > N/2, in other- words, the requirement that the antisymmetric 
component of the spin representation is more than half-filled. 



C. Calculation of the Magnetization using Gaussian Fluctuations 

One way to examine the consequences of this residual coupling on the single-ion Kondo 
effect is to compute the field-dependent magnetization of the ground-state. The application 
of a field provides a controlled way of examining the cross-overs associated with screening 
processes. To compute the magnetization, we need to introduce a magnetic field and calcu- 
late the field dependent ground-state energy, including the effect of the Gaussian fluctuations 
around the mean- field theory. For SU (N), there are N — 1 ways of introducing the magnetic 
field. We shall use the form 

H B = -B^m^infa + n bc ) (63) 

a 

where m\ = —1, m 2 = 1. With this choice, the field splits off two bosonic and two fermionic 
states from the other (N — 2) levels. 

The mean-field free energy in a field is then given by 

NV 2 

+ —f-\Qo-(Y (64) 



F mft [B] = Y: 



where £ CT = £ — Bm a , A = A — Bm a To calculate the magnetization, we must differentiate 
the Free energy with respect to B. ( Since the Free energy is stationary with respect to 

19 



changes in A and (, we do not have to worry about how these fields change w.r.t. B). The 
magnetization is then 



M = n bi + (n/i - ra/2 ) 

m f (B) 



2S+- 



-if A \ -i ' A 

tan — — — tan 



(65) 



— -B — A/ J \B — A/ / 

The first term is the residual unquenched moment, the second term represents the spin- 
polarization of the Kondo singlet. Technically, we should lump this term with the other 
0(1) corrections that we need to calculate from the Gaussian fluctuations. The mean- field 
only reliable predicts the terms of order N, and thus to this order, 

M = 25 + 0(1) (66) 

To calculate the 0(1) term, we need to include the zero-point corrections to the ground-state 
energy. 

There are two types of Gaussian fluctuation around mean-field theory: bosonic fluctua- 
tions in V, A and (, plus the fluctuations of the a field. Fluctuations in V and A/ = A + ( 
are associated with the interactions in the Fermi liquid. These terms renormalize m,f(B) 
and produce an order 0(1/ N) correction to the magnetization. Fluctuations in A& = A — ( 
renormalize the entropy of the free moment, and do not produce any correction to the un- 
screened moment. The only 0(1) corrections to the magnetization are those associated with 
the fermionic fluctuations. The Lagrangian for these fluctuations is given by 

Hi 

C a = Q a\d T + 2()a + \f%(d T + 2()a + a\d r + 20&U] (67) 

In a field the number one boson condenses, so that Af,i = 0. To fulfill the constraint on the 
Fermi fields, n*j = constant, we require that A/ = A&i + 2£ + B is field independent, which 
implies that in a field, 2£ = 2( a — B. 

When we expand the effective action S e ff to Gaussian order in the afield, we obtain 
the correction 
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where 



S a = - J dld2a{l)D-\\ - 2)a t (2) 



D-\l - 2) = Q D (d T + 2C) - (d T + 2C) 2 (T^(1)^(2)) 

= Q (d T + 2C) - (d T + 2C) 2 £ G /tr (l - 2)G 6ct (2 - 1) 



(68) 



(69) 



is the inverse propagator for the a fermion. Written out both diagrammatically, and in the 
frequency domain, this is 



(70) 



D~\u) = Q {u-2Q - N[t(2( - lu) 2 ]<S>(u) 



where 



T 



It proves convenient to factorize D 1 (cu) as follows 



D-\uj) =Q (uj-2QP{uj) 



(71) 



(72) 



where 



P(uj) = 1 + ^—K^) 



A detailed calculation of P{oj) at zero temperature ( Appendix D) yields 



P(uo - id) = P {uo + B-iS) 



(73) 



(74) 



where 



P (u) = A f (u) 



u-2Co L T K 

A In h nxrib 

ixq \ uj 



(75) 



is the zero-field expression for P{uj) and Af(u) = ImGf{uj — iS) is the zero-field f-spectral 
function. 
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Now the free energy associated with the Gaussian fluctuations in a is given by 
F « = -Tj2{H-D-\itu n )] - ln[Q Q (2C - iu n )]} e^ 0+ 

n 

= -T^ln[P(^ n )e-» 0+ ] 



(76) 



where we subtract the Free energy of the auxiliary a field in the absence of interactions, 
to avoid over-counting. Carrying this sum out by contour integration, and taking the zero- 
temperature limit, the zero-point energy is 



/U H U 
— lmln[P (uj + B - iS)} 
-D 7T 



(77) 



where we have inserted the field dependence by replacing P(u>) — > P {w + B). Differen- 
tiating this with respect to the applied magnetic field, the screening contribution to the 
magnetization due to the interaction between the spin and Fermi fluid is given by 

dB 
r° du 



— dJm\n[P (u + B - i5)} 

-D IT 



1 i° 

-lm\n[P (u + B -iS)] 

7T 



D 



(78) 



The lower limit of this sum gives a contribution — ^Imln[— D \n(D/T K ) + n b )] = —1. 

The final result for the fluctuation contribution to the magnetization is then 



M a (B) = -1 + -1mm 

7T 



(79) 



ixq \ B 

To obtain the total magnetization, we must add this to the result M m f t = 2S + rrif(B) 
obtained from the mean-field theory. The total magnetization, evaluated to order 0(1) is 
then 

Mi(B) M 2 (B) 



M(B) = [2S - 1 + m f {B)} + -1mm 

7T 



nq 



B 



(80) 



where 



1 f 


"2C + B' 




\2(o-B\-\ 


— < arctan 


— arctan 


7T { 


A 




[ A iJ 



(81) 
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FIG. 4. Showing the field dependent magnetization for a Kondo model with ( > 0. In this 
case, the Kondo effect is a single-stage process. Inset shows the corrections to the magnetization 
derived from the residual interaction between the heavy fermi and spin fluid. These corrections 
are positive (ferromagnetic) at low temperatures and remain small at al temperatures 

We can check this result by completely removing the fermionic contribution. In the limit 
A — > 0, 2( Q — > T K , Mi(B) and M 2 (B) develop discontinuities which precisely cancel, to 
yield 



M(B) = 2S --{ arctan 

7T [ 



Ik 
B 



2Sn 



7T 

+ 2 



(82) 



which is the residual magnetization of the spin S Kondo model. 

Let us now restore a finite A. The first term, Mi(B) represents the screening of the local 
moment, by the Kondo effect. Mi(B) has the limiting values 



M l (B) 



2S-1 (B«T K ) 
2S (B » T K ) 



(83) 
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FIG. 5. Field dependent magnetization for a Kondo model with £ < 0. In this case, the 
Kondo effect is a two-stage process. The corrections to the magnetization derived from the residual 
antiferromagnetic interaction between the heavy fermion and spin fluid. These interactions grow 
to —1 at low temperatures, corresponding to a second screening process at the renormalized Kondo 
temperature T^. 

The second term M 2 (B) derives from the screening effect produced by the residual interaction 
between the Fermi liquid and the magnetic moment. For ( > 0, this can be re-written 



M 2 {B) = -arctan 

7T 



2Sn 



Ik 
B 



(84) 



At low fields, this contribution is small and positive, corresponding to an irrelevant residual 
ferromagnetic interaction. At 2( = B, this contribution passes continuously through zero, 
due to a change of sign in the residual interaction. At still higher fields, this correction 
remains small, and asymptotes to zero. 

By contrast, if ( Q < B, in this case the residual interaction with the Fermi fluid is 
antiferromagnetic. Since 2( Q — B is always negative, we can now write 
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For small fields, M 2 (B) — ► —1, so that this term constitutes an additional screening contri- 
bution to the magnetization. For small negative zeta, we can approximate this expression 
by 



where 



M 2 {B) 



1 



arctan 



7T 



+AHn 



B 



2nS 



(86) 



T,* rji — Z!* | rji — 7T(7COt[7rfn^— h) 



(87) 



which corresponds to a second screening process, governed by the second-stage Kondo tem- 
perature T^. 

An alternative way to derive the same result is to consider how the Schwinger boson field 
condenses in an applied field. The constraint associated with the bosonic part of the spin is 
written fllHD 

f 



2S = n b + 



Q 



(88) 



In zero field, the Fermi fluid is unpolarized, and the magnetization is given by the condensed 
part of the Schwinger Bose field. Suppose we apply a small field that condenses the b\ 
component, so that 

(89) 



K = VMS al + Sb 



then the constraint can be re-written as 



2S = M + y £(5blSb tr ) + - 

a V 



(90) 



or 



Q 



(91) 
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Diagrammatically, these two contributions to moment reduction are given by 




Of course, only the combination of the two terms is gauge invariant, but by fixing the gauge, 
we can assign them each physical meaning. The first corresponds to fluctuations in the 
direction of the local moment. The second represents the reduction in the amplitude of the 
moment derived from the inter-conversion of spins into heavy fermions. We can compute 
the sum of these diagrams by noting that they are generated by differentiating the RPA 
diagrams contributing to the fermionic zero-point energy with respect to the frequency: 




which enables us to identify the reduction in the magnetization with the fluctuations in 
direction and magnitude of the local moment. In this way, we see that the fluctuations 
which screen the moment are given by 

^ 1 [ 1, (n* f < N/2) 

^ I 2, (n* f > N/2) 

depending on whether a one, or two-stage screening process takes place. Although these re- 
sults are only calculated to leading order in the large- N expansion, we expect the appearance 
of integer values for the screening is exact for a local moment. 
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IV. STRONG-COUPLING PICTURE OF THE TWO-STAGE KONDO EFFECT 



To gain a complimentary insight into the two-stage Kondo effect, it is useful to examine 
this phenomenon in the strong-coupling limit. Imagine a local moment, described by an 
L-shaped representation of SU(N), denoted by the Young Tableau 

— 25 — 



(96) 



In the ground-state of the strong-coupling Hamiltonian 



(97) 



electrons form a singlet with the fermionic part of the spin creating a partially screened 
moment, denoted by a Young-tableau with a completely filled first row. 

- — 25 — 

25-1 



s* = (r eo + s) = 



N 



(98) 



where in this example we have taken N = 8. Since the first column of the tableau is a singlet 
(with N boxes), it can be removed from the tableau, leaving behind a partially screened spin 
5 — 1/2, described by a row Young-tableau with 25 — 1 boxes. If we now couple the 
electron at the origin with electrons at site T' via a small hopping matrix element t « J, 
then the virtual charge fluctuations of electrons in and out of the singlet at the origin will 
lead to a residual coupling between the partially screened moment and the electrons at the 
neighboring site '1' 
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J* 

H(\) = J-jr S* ■ C la r a pC^ 



(99) 



where J* ~ t 2 / J. In the SU(2) Kondo model, only electrons parallel to the residual moment 
S* can hop onto the origin, which gives rise to a ferromagnetic coupling J* < 0. In the 
SU(N) case, electrons can hop provided they are not in the same spin state as electrons 
at the origin. The sign of the coupling J* depends on the number of conduction electrons 
n c = n — n*j-, bound at the origin. If the n c = N — 1, electrons hopping onto the origin will 
have to be parallel to the residual spin, so in this case the coupling is ferromagnetic, J* < 0. 
By contrast, if n c << N, there are many ways for the electron to hop onto the origin with 
a spin component that is different to the residual moment, so the residual interaction will 
be antiferromagnetic, J* > 0. By carrying out a large N calculation in the strong coupling 
limit or by making a detailed strong coupling calculation for SU(N), we are able to confirm 
that for N > 2, J* changes sign when the number of bound-conduction electrons is less than 
N/2, and in the large N limit, it is given by@ 



J* 



t 1 



J{l-n))n) 



n 



f 



n) + n b ) 



(100) 



where n* f = n* f /N and n b = 2S/N. 

When n*j > N/2, 4 J* > 0, the strong-coupling fixed point becomes unstable, and a 
second-stage Kondo effect occurs, binding a further N — 1 electrons at site "1" to form a 
state denoted by the tableau 

- — 2S^ 



(r eo + r ei + s 



N 













Cl 






Cl 




Cl 




Cl 


c c 


Cl 


c c 


Cl 
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Cl 



2S-2 



(101) 



corresponding to a residual spin S** = S — 1, M = 2S — 2. This final configuration is stable, 
because an electron at site "2" can only hop onto site "1" if it is parallel to the unquenched 
moment. 
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We see that our supersymmetric approach permits us to examine the consequences of 
this two-stage Kondo effect, starting from weak coupling. Translated into the weak coupling 
language, the two vertical columns of the Young tableau will correspond to two separate 
screening clouds, of very different radii 

l = £, and r=£, (102) 

respectively. It is remarkable that a point-like complex impurity can give rise to two separate 
length-scales in this way. 



V. DISCUSSION 

In this paper we have developed a spin representation that interpolates between the 
Schwinger boson and Abrikosov pseudo-fermion representations, which exhibits the prop- 
erty of supersymmetry. As an exploratory exercise, we have applied the method to to a class 
of single impurity Kondo models, where we have been able to examine how local moment 
behavior emerges as the strength of the Hund's interactions between the spins is systemati- 
cally increased. Suppose we consider a spin representation with Q boxes, and examine how 
the ground-state evolves as we progressively increase the moment from S —1/2 to S — Q/2. 
One of the surprising discoveries, is that there are in fact, two routes by which the magnetic 
moment emerges from the Fermi liquid (Fig. 5): 

1. One stage Kondo effect, where Q < N/2. Once the spin S exceeds one half, a partially 
screened moment is generated. The low temperature fixed point is described by a co- 
existence of a Fermi liquid and a moment of spin S* = S — 1/2, with a slowly vanishing 
ferromagnetic coupling between the two degrees of freedom. 

2. Two stage Kondo effect, where Q > N/2. At intermediate scales, the moment quenches 
to a spin S* = S — 1/2, but the residual coupling to the conduction sea is now 
antiferromagnetic, and a second-state quenching occurs to a spin S** — S — 1. When 
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(a) Q< N/2 One stage Kondo 



FL 



2S 1 



Under- screened S =S-l/2 
2 3 Q-2 Q-l Q 



(b) Q>N/2 



2S 1 



S = S-1 



Two-stage Kondo 



S = S- 1/2 
Q-2 Q-l Q 

One stage Kondo 



FL FL 

FIG. 6. Two scenarios for the emergence of magnetism as the size of S is progressively 
increased, (a) Single Stage Kondo effect, where for S > |, a partially quenched moment with 
spin S* = S = \. (b)For Q > N/2, the initial emergence of a local moment is accompanied by a 
two-stage Kondo effect, where, provided n*c > N/2, the spin is screened from S to S** = S — 1. 

the starting spin is S — 1, a new singlet phase is formed, with one additional fermionic 
bound-state. ( We label this phase FL* in Fig. 5, but do not at this stage know if 
this state is a Fermi liquid.) 

We should now like to discuss the future extension of this approach to a Kondo lattice. 
Two decades ago, DoniachB argued that the properties of a Kondo lattice should depend 
critically on the ratio of the Kondo temperature to the RKKY interaction k = Tk / Jrkky- 
Heuristically, the Kondo and RKKY scales are related to the Kondo coupling constant 
according to 



T K ~ De-V Jp , Jrkky ~ J 2 / D 



(103) 



where D is the conduction band-width. Doniach pointed out that k grows with the size 
of J, arguing that for small J, the system is expected to be antiferromagnetically ordered, 
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but for large J, the magnetism melts to form a "heavy" Fermi liquid. Unfortunately, our 
theoretical understanding is at present, only limited to a discussion of energy scales, and 
little is known about the nature of the transition between these two limiting cases. 

Can we do to shed light on these issues by extending the supersymmetric approach to 
the lattice? The mean-field solution will in general give rise to a heavy Fermi liquid which 
co-exists with a lattice of under-screened moments. The Fermi surface volume Vfs will 
depend on the size M of the under-screened moments, 

^ = (n c + n))/N = (n c + Q- M)/N (104) 

At the mean-field level, these moments can point in any direction, but once we include the 
effects of the Gaussian fluctuations of the a fermions, two effects will take place: 

• The local moments will be partially screened by the Kondo effect with the heavy 
electron fluid. 

• The fluctuation free energy will be sensitive to the direction in which the spins con- 
dense. 

Typically, we expect the fluctuation free energy will be lowest in an antiferromagnetic spin 
configuration, where, for instance 



(6 ff (x)) = V2M[cos 2 (Q • x/2)5 CTl + sin 2 (Q • x/2)5 ct2 ] (105) 

This dependence on the relative orientation of the moments defines a renormalized "RKKY" 
interaction J* RKKY - By tuning J we will be able to examine how the RKKY interaction is 
renormalized by the Kondo effect and examine how the staggered magnetization depends 
on the screening process 

M = 2S-m(J) 
m(J)=T,(^Jb a ) + ^9) (106) 

cr V 

where m(J) is a continuous function of J. The critical value where M vanishes (at small 
S) defines the point where the magnetism is eradicated by the Kondo effect. In this way, 
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Y=0 



Q =1 [? 



Y= -4 




Y= -2 



Y=0 



Y=2 



Y=4 



FIG. 7. Illustrating the idea that the properties of a class of L shaped representations will 
enable us to triangulate on the properties of a small N Kondo lattice. 

we hope to be able to carry a model calculation of the phase diagram first envisaged by 



An open question, is whether this new approach can shed light on the nature of the 
quantum critical point separating the magnetic and the paramagnetic phases of the Kondo 
lattice. Ultimately we are interested in the properties of a Kondo lattice of elementary 
spins corresponding to one box in Young Tableau. This large N approach can provide 
information about the properties of a class of L-shaped representations. If there is any 
universality associated with the emergence of magnetism at absolute zero, then perhaps a 
large N approach will enable us to triangulate on the properties of a real Kondo lattice. 



One of the interesting aspects of the supersymmetric approach, is the appearance of 
fermionic "phase fluctuations" between the spin and the heavy electron fluid: it is these 
fluctuations, described by the gauge fermion a, which mediate the interaction between the 
magnetic condensate and the heavy electron fluid. In the case where the "mass" of this 
excitation, ( is positive, the gauge fermion a can be integrated out of the problem, and the 
interaction between the magnetic and electron fluid could be treated as a point interaction. 



Doniach. 



(Fig. 6) 
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However, when ( is negative, the a fermion gives rise to a new bound-state. Will the same 
phenomenon occur in the lattice? This could lead to the possibility of two different kinds of 
fixed point: 

• A Millis Hertz fixed point!!, where the weak ferromagnetic interaction between the 
magnetic and electron fluid can be treated as a point vertex. In this situation, the 
transition will be described by the interaction between a Gaussian magnetic fluid and 
a well-defined Fermi surface. 

• A non- Fermi liquid fixed point, where the dynamical fermion mediating the magnetic 
interaction becomes an active participant in the physics. If the gauge fermion de- 
veloped gapless excitations, then the decay of heavy fermions into unquenched spins, 
described by the process 

+ (107) 

would lead to a phase with a novel kind of spin-charge separation. 

These points will be examined in greater detail in a future publication. 

We should like to thank Natan Andrei, Antoine Georges, Mireille Lavagna, Olivier Par- 
collet and Revaz Ramazashvili for critical discussions related to this work. This research 
was supported in part by NSF grant DMR 961999 and research funds from the EPSRC, UK. 

Appendix A: Operator Expression for the Cazimir 

In this section, we prove that the Cazimir 

S 2 = (tfVb + frf) 2 (108) 

can be written 

S 2 = Q(N-y-jt). (109) 

To show this relationship, we use the completeness result. Using the normalization 

33 



Tr[FT 6 ] = 5 ab (110) 



this is 



**** + ±5 ab 5 cd = 5 ad S bc (111) 



By using this to expand 



S 2 = (tfTb + pTf) 2 (112) 



we obtain 



S 2 = -JjQ 2 + (b% + flhMK + flfa) (H3) 
We can expand each term in this expansion as follows: 

b%b\b a + flf b flf a = n b {n b + N-l)-n f (n f -(N+ 1)) 

= n 2 b -n) + NQ + (n f -n b ) (114) 

Also 

flf b blba = -blf b fX + n b (115) 

Combining these results we obtain 

S 2 = -^Q 2 + (n 2 b - n)) + (N + 1)Q - 2b\hf a b a 
= -^Q 2 + (N + l)Q + Q[n b -n f - ^} 
= -^Q 2 + NQ- Q[n f -n b + 6*]] 

= Q(N-y-jj) (H6) 

where we have used the result 299^ = [9, 9^] + Q to carry out the last step but one. 
Appendix B: Supersymmetry of Lagrangian 

In this section, we examine the transformation of the Lagrangian 
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under the transformation ~^> a — > h^> a , where 

h = e 1 ^ 



C 0a 
'at 



V 1 _ W »7 



e *e c q 
e" ie c 



-7] Vi - 7777 

is a general member of the group S"?7(l|l). Under this transformation, 

C 0<*" 



C^-C^ + E^U^ + ^t 

(T » 

When we expand the correction, we obtain 



)at -C 



E*U^)* CT = E*t 



%Qq + i0 f T 3 + T]d T f] + 



e 2ie id T r] 



—e c d T r] 



2Cr/e 



and 



(e 2id ( - l)0oJ fj0a + {0a*)ri 



Combining (|i~20D , and (|i~22l) we obtain 



£:«*-£*£ ^ + A' + 



117) 



;ii8) 



(119) 



(120) 



'12U 



(122) 



(123) 



where A' = A + z'0q, £' = C + > a ' = e_;i6,c ( a + ^) ; a '^ = e l9( (a^ + ff). The primes on the 
gauged derivatives denote = (d + 2(') and $ = (—9 + 2('). The Lagrangian is thus gauge 
invariant under the transformation 

A — > A — idq, (^(-ie 



(124) 
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Appendix C: Evaluation of fermionic Mean-Field Free energy 

We wish to calculate 

F f = —NT^\n[\f + iA n — iuo n ]e iuln0+ (125) 

n 

where A n = — Asignw n . We shall regulate this sum by calculating 

F f = -NT*£(\n[\ f + iA n - iu n \ - [X f -> X f + D])e l ^ 0+ (126) 
which we re-write as 

a0+ 



Ff = -2NTReJ2(Ht + iuj n ] - Ht + D + ^n])^ 



n>0 



= -2iVTRe^(ln 



£ 1 
— + - + n 



n >^ ^iT 2 
where £ = A/ + iA. Next, using the result 



In 



Z + D 1 

— H h n 

2mT 2 



e m0+ (127) 



n>0 

this becomes 



£(ln[& + n] - ln[a + n])e m0+ = In (^j) + - a) (128) 



^ = — 2TRe In [2 7 T 2 - - (129) 



Appendix D: Calculation of P(uo) 
We begin by writing 



where 



P(u) = 1 + [ ^—K) $(o;) (130) 



# , 

To calculate $(u;), we replace the discrete Matsubara sum by a Contour integral, to obtain 

$H = ^ E / ^n(z)G bir (z)G fa (z + W ) (132) 

36 



where the integral runs counter-clockwise around the poles in the Green's functions. Using 
the spectral decomposition, 

G fa (z) = f ^A fa (e)-±- (133) 
J 7T z — e 



this becomes 



*M = ~ E / -K + m\ — l —^- A f^) (134) 



where = n(A& CT ) is the Bose occupancy. Now 
/ - -—T XT A /-( e ) = E 



2 9 

^-2 C + 8 A„ =2SG ^ + i;) <135) 
where we have replaced 2( — > 2( Q — B, and Gf(icu n ) = (iu n — 2( + iA n )~ 1 is the f-propagator 
in the absence of a field, so that 

P(u) = 1 - ^ + ^~ 2C ° [G f (u + 5)n 6 + /] (136) 

where n b = 2S/N and 

where -A/ (a;) = ImGf(uo — iS). Now since AT — 2 of the levels are unshifted, to leading order 
in the large iV expansion, we can set m a = in this expression. Also, since A& = B in a 
magnetic field, we can write I = I(uo + B) where, at T = 0, 

/(„) = r &mo j_ (138) 

J-d 7r a; — e 
Combining these results together, we can write 

P(u, B) = P (lo + B) (139) 

where 
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P (co) = 1 - [G f (u)n b + I(co)\ 



(140) 



is the zero field form of P(tu). Going on to evaluate I (to), we obtain 

e 

1 



r° A f (( 

J-D 7T 

f° *L ( 1 

7-u 27ri \e — 



c<j — e 



1 



rfe 



1 1 \ 1 



d 2ni I \ e — £ e — to I £ — to 



27ri I — £ 



1 .„ i 



1 In (£ 



a; 



a; 



(141) 



so that 



1 ~ 2Co) 
9 



n fe 1 



1 .„fi 



(142) 



Now by writing £ = T K e mnf = £* e ' t27rn f , where n/ = n*f/N, we can put this in the form 



PoM = i 



(w - 2C) 



^ + ii„dU l_ 



Since rib + rif = q, there are cancellations between the first two terms which give 



u-i) nq 
Another useful way to rewrite this expression, is 



tO 



(143) 



(144) 



Po(w) 



-iA 



A / 1 
—Re 

qn 



In £ 



tO 



(145) 



To make contact with the bosonic Kondo model, it is useful to split the first term into a 
real and an imaginary part, so that 



-iA 



AA f (u) - iARe 



One can then move the second term above into the logarithm, writing 



-iA 



Re 



qir 
to obtain 



In (^-\ -iARe 



-iA 
qn 



Re 



1 ^ K -L ■ ~ 

in h 27m b 
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P (u; - i<5) = 



to -2( f T K 
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ln h i7rnb 



(146) 



(147) 



(148) 
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